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ABSTPvACT 

An   analysis   is   done   of  the   motion   of  an   underwater  cradle    for  submerged 
recovery  of  the   DSRV  Alvin. 

The   cradle   is   modeled  as    a  three   degree   of   freedom  system  with   four  cor- 
ner springs   and  Includes  heave   pitch   and  roll  motions.      The   surface   ship   Lulu 
is  taken  as   the   driving  force   that   couples   energy   into  the    cradle  through   the 
springs.      Hydrodynamic  damping  values    are    found  by  a  54  section    finite  ele- 
ment  analysis   and  are    compared   to   corresponding  linear  approximations. 

Solution   of  the  resulting  equations  was   accomplished   using   a  Runge-Kutta 
numerical  Integration  m.ethcd   developed  at  M.I.T. 
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Ibf  Will  substantially  decrease  the  transmit ^ed  ship  motion.   Adaitionaiiy , 
the  results  indicate  that  the  pitch  and  roll  natural  frequency  is  higher  than 
the  heave  natural  frequency  and  consequently  is  the  limiting  design  considera- 
tion. 

The  results  of  the  analysis  are  discussed  and  several  possible  approaches 
are  considered  for  producing  a  soft  spring  constant  while  at  the  same  time 
supporting  the  high  static  load  of  the  cradle  weight. 
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CHAPTER   I. 
INTRODUCTIO:f 

The  handling  of  heavy  objects   through   the   air-sea  interface   is  not   a 
new  or  unique  problem  and  has  been   studied  many  times   in   the   past.      The 
deep   s submergence   research   vehicle    (DSRV)    Alvin    owned  by   the   U.S.    Navy   and 
operated  by   Woods   Hole    Oceanographic   Institute    (V/HOI)    is   most    likely   the 
DSRV  with   the  world's   best  record  for  successful   dives    completed.      The 
submersible   is  handled  by  its   support  ship,    a  catamaran  named   Lulu  and  is 
transported  to  and  from  the   open   ocean  dive   site  betv;een   Lulu's   twin  hulls. 

The   launching  and  retrieving  operation   is   accom.plished  with   a  great 
deal  of  care,   Alvin  being  initially  positioned  on  her   30-ton   capacity   cradle 
just   aft   of  amidships.      The   cradle,    (Figure    1),   is   supported   at  its    four 
comers  by   chain  hoists  which   are   two-blocked  such   that  the    cradle   is   at 
deck   level  during  transit.      This   allows   easy  maintenance,   repairs   and  m.an- 
ning  of  Alvin. 

During  launching,   the   catamaran   lies  with   its   stem   to  the   seas,   the 
cradle   is   then   lowered  and  the   floating  submersible    clears   Lulu  with   the   aid 
of  Alvin 's   screw  and  numerous   snubber  lines.      Recovery   is   the   reverse   of  the 
operation,    Alvin  surfacing  well  away   from  Lulu,    closing   Lulu,   and  finally 
driving  between   Lulu's  hulls   and  being  manually  positioned  over  the   sub- 
merged cradle.      Last,  when   the  wave   timing  is    correct,   the    cradle- 
siibiiiersible    combination   is   raised   clear  of  the  water. 

This  method  of  launching  and  retrieval  has   allowed  reasonably  routine 
operations   to  occur  in   sea  states    up   to  three,  with   occasional  sea  state 
four  retrieval  due   to  emergency   conditions. 
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Figure   1 . 
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V/it!;out   going   into   a  discussion   of  sea   state   probabilities,    it   is   v;ell 
knov;n   that   along  the   north   Atlantic   off   Cape    Cod   the   sea  state   v;ill   frequent- 
ly  go   fron  state   three   to  state    four  due    to  a   prevailing   southv/est  wind  which 
increases    fron  ten  knots   in   the  morning  to  about   20   in   the   afternoon.      Sea 
state    four  makes   for  a  very   difficult   recovery  should  it   occur   during  the 
time   of  a   dive.      The   problem  here    really   is    that   once   the   Alvin   is    in   posi- 
tion  to  be   hoisted  aboard,   the    relative   motion  betvieen   Alvin,    the   cradle,    and 
Lulu  could   cause   very  high   and  damaging  shock    loading  on   Alvin,   the  hoist, 
or  the    cradle.      This    loading   could  possibly   result   in  severe   damage   to  Alvi© , 
or  even   another   loss. 

Clearly  then ,   a  smooth    docking  to  a  stationary   cradle    followed  by   a 
rapid  hoisting  through   the   air-sea  interface  would  be  highly   desirable   and 
wniilfi   vi-ptiiailv  eliminate   shock   loading  and   consequent   damage. 

Fortunately,   the  Alvin 's   cradle   consists   of   four   120'    lengths   of  chain 
and   four  synchronized  wildcats.      If  the   cradle  were   lovjered  to  approximately 
100'   belov;  the   surface,   the   surface  wave   effect  would  be   virtually  inefec- 
tive   on   Alvin  and  it  would  be  possible   to  dock   or  depart   from  the   cradle 
without   the    fear  of  poor  surface    conditions.      The  problem  which   appears   to 
be   the   limiting  factor  is   the  motion   of  the   cradle   itself.      This  motion   is 
nearly  the   same   as   that   of  Lulu   due   to  the   inelastic   chains  which   connect 
the  two.      Some  work  has  been  done   to  show  the   transmitted  surface  motion 
that  the   cradle  would   feel.      It  has  been   generally   concluded  that   if  members 
more  elastic  than   chains  were   used,   perhaps   the    cradle   motion    could  be  un- 
coupled  from  the   Lulu  and  the   Alvin   could   dock  v;ith   a  virtually  motionless 
cradle   at   100'    and  then  be  hauled  to  the   surface   rapidly   and  safely   through 
the   air-sea  interface. 
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T!io   complete   analysis   of  the    Lulu-cradle    interaction   is    the   priir;ary 
;oncern   in   this   thesis. 
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CliAPTSK  II 


T!;eopy  a::d  devclopment  of  eouatio;;s  of  motion 

To  this   tine   the   ship   cradle   interaction   rr;otion  has   been   studied   using 
a  simple   model  shov;n   in   Figure    2.      While   this   model   applies    very  well   to  a 
svstem  v.'ith   a  single   support   it   does   not  prerlict   the   motion   of  a  four  point 
suspended   cradle  with   ship  pitch,   roll,   and  heave. 


SHIP 


X^=Yosir.wt 


K 


M 


Ity 


77777777777777777-17 


y(tj=yQSinut 


Figure  2 . 


Looking  at  Figure  3  it  is  apparent  that  the  motions  of  a  cradle  be- 
neath a  ship  are  extremely  complex.   The  cradle  can  exhibit  vertical  motion 
(y),  pitch  (9),  roll  (6),  fore  and  aft  sway  (a),  transverse  sway  (6),  tv:ist- 
ing  ibout  the  y  axis,  as  v:ell  as  horizontal  motion  in  the  x  direction.   In 
efrect,  all  six  degrees  of  freedom  can  be  achieved  v:ith  the  cradle.   For- 
tunately, to  soxewhat  simplify  the  problem,  fore  and  aft  sway  (a). 
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Figure  3 
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trannverse   sv/ay   (B),    ano    tv;ist  were   oxperimental.ly    found  to  be   insign5.fi- 
cant. 

Consequently,   the  problem  to  be   solved  is   a  three   degree   of  freedom 
one  v/ith   roll   (<>),   pitch    (9),    and  heave    (y),   and  where   all  three  modes   of 
motion   can  transfer  energy  back   and   forth   thus    contributing  to  the   com- 
plexity  of  the   problem. 

Before  beginning  the   modeling  the   initial   conditions   must  be   defined. 
To  begin  with,   the   cradle   is   assumed  to  be   of  length   1  =    17.5    feet,  width 
w   =    15.3   feet,  with   a  v/eight   of  approximately    508   slugs.      Since   the 
reference  material  has   indicated  that   a  possible  mode   of  operation  would 
be  with   a  flat  plate   affixed  to  the  bottom  of  the    cradle ,   this   is   the  mode 
chosen   for  our  model.      Consequently   the   cradle    is   assumed  to  be   a  rectan- 
gular  flat  p]ate  with   ma-^s   equivalent  to  that   of  th^    cradle   alone.      In   the 
final   analysis   it  will  be   shovm   that  the   added  mass   terms   associated  with 
the   cradle   so  overpower  the   actual  mass   that  the   exact   configuration  of 
the   cradle  structure   is   inconsequential. 

It  was   first  necessary   to  pick   an   adequate    co-ordinate   system  which 
would  accurately   describe   the  system  and  make   it  possible   to  completely 
position   the   cradle  knowing  only  three   co-ordinates — heave    (y),   pitch    (e), 
and  roll   ((^).      A  reference   system   (Figure   U),   fixed  above   the   cradle  was 
chosen  with  X  positive    for\-;ard,  Y  positive   down   and  Z  positive   to  port. 
Using  the   conventional  right  hand  rule  ,  positive   9    and  positive  ^   are   as 
shown   in   Figure   4.      As   is   conventional   in   dynamic  analysis   the   cradle   is 
given  an  instantaneous   displacement  positive  y,   positive   6,   and  positive   (^. 

Since   the  equations   of  the  system  which   are  being  dealt  with   are 
quite   complex  vvith  many  terms,   it  will  be   convenient  to  set   up   the 
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CradlG  Equilibriim 


Figure  H 
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resultant  equation   in   standard   state    variable    form   c^onsisting   of  six   first 
order  dif fer'^nntial  equations    and  an   output  expression.      The  most  general' 
form  of  such   an  equation    is 


Xi(t) 
X2(t) 


X„(t) 


All  Ai2 hn        ^l(^) 

A21   A22 X2(t) 


Anl  An2 ^^n        X^Ct) 


b lib 12 


bnl- 


••bin  Ul(t) 


.bnr         U„(t) 


or 


X(t)   =   AX(t)    +   Bu(t) 


where   A  is   the   N   x  N  system  matrix  with   B  an  N   x   r  control  matrix. 
btrictiy  speaking,   tnis    above  representatioii   applies    only   to  a  N ""'  oraei' 
linear  equations ,  however,   since   there  will  be   a  multitude   of  terms,  both 
linear  and  nonlinear,   it  will  be   a  good  bookkeeping  aid  to  arrange   terms 
as  such.  .        • 

The  basic  equation  of  motion  for  a  spring,  mass,  and  damper  system 
is  used  here.  Development  of  the  equations  of  motion  is  begun  with  the 
consideration   of  this    free  body   diagram. 
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I   forces   =   ma 


V7EIGHT 


^y  -  "  ^drag  ~  ^  spring  +  force  (exciting) 


and  is  the  case  of  three  degrees  of  freedom 
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^2'^  ~        ^drag   ^spring  '^   ^exciting  ("^o^ 
Each  term  will  be  discussed  independently. 

The  Soring  Force 
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Figure   5 


The   four  springs    (chains)    at  the   comers   of  the   cradle   are   the  major 
factor  in   decoupling  the   surface   motion   from  the   cradle.      The   spring   force 
is  proportional  to  its   stretch   and  this   stretch  may  well  be   different   for 
each   of  the    four  springs.      This   is   one   of  the  mechanisms  by  which  energy 
is  transferred   from  one   mode,   such   as  heave,   to  pitch   or  roll.      Ideally, 
if  all   four  springs  were   equal  and  the   cradle  were   given   a  pure   pitch,   roll. 


or  heave  motion,   then  the   other  modes  v;ould  not  be   excited  and  you  would 
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havG   s.-iir:p..ly  throo    uncouoled  r:ystcnis.      However,    considering   the    inexactncoS 
of  thf:   Real  V/orld,   the    croGs    coupled  damping,   and   the  exterior   forces,    it 
would  be    unrealistic   to  perform  such   a  simple    analysis.      In    actuality   the 
ship  pitches,   rolls,   and  heaves   about   an   equilibrium  position   as   does   the 
cradle.      The   stretch   of  the    four  springs    is    defined   as   the   difference   in 
position  betv.'een   the   position   of  the   ship   and  the   position   of  the    cradle 
relative    to  their  respective   equilibrium   for  each   motion   mode. 

The   respective   stretches   of  the    four  springs   K,  ,   K2 ,   K3 ,    and  K^   are 
as   follows :  ■  • 

^l  =  (y-YQsinu;Lt)  +  (L/26-L/2eQsina3j^t)  +  (w/24'-w/2$Qsina32t ) 

.62   =  (y-YQsina)j_t)  +  (L/28-L/20Qsinw-Lt)  -  Cw/2(j)-w/2<l'Qsinw2"t) 

63  =  (y-rQsmcOj^t)  -  (i,/2^-L/2^QSxnuiit)  -  (W/2<p-w/2$QSinw2t; 

64  =  (y-YgSinw^t)  -  (L/2e-L/2e   sinw-Lt)  +  fW/2({)-w/2$sinu)2t) 

This   is   assuming  thatYo,e„,   and  $„ ,   are   the  ship's   pitch,   roll,   and 
heave . 

Carrying  this    a  bit   further,  when   the   stretches   are  multiplied  by 
their  spring   constants ,  we   get  the   respective   restoring  spring   forces   and 
moments   in   general   form. 

The  equation   for  heave — 

my  =   -   K^6^  -  K^5^   -   K363   -   K^\ 
The  equation   for  pitch--  ' 

1^6   =    -   K^l^^   -    K2|<S2    +   K3|^3   +   K^^^^ 
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The   equation    for  roll-- 


W  w  W  w 


or  plugging  in    6's, 
heave  equation-- 

.     my   -h  Kj^Sj^  +   K262   -i-  K363  -1-   K45^   =   0 

Plugging   in: 

^^2   '   ~   ^I'-^^l   ~   YQsinw^L"^)    '''    ^o'^S  ~    2    Qo^^^^l"*^^    "*"    ^^  ^5    "   2"  "^0   3inw2t)] 

L  L  WW 

-  K2[Cy2   -   YQsinco-|_t)    +    ^o"  ^3   "   "  Qq^^^'^I''-^    ~    ^~  ^5   "   2"  ^q^'^^'^S"'^^^ 

L  L  WW 

-  l^sC^Yl   -   Y^sinco^t)    -    i^^Z 8   sinw-^t)    '    ^2^^    ~  'l    ^q^^'^^2^^'^ 

pitch  equation-- 

L  L  L  WW 

I^yil  =   +  K3  ^(y-L   -   Y^sinu-Lt)    -    (^  Ya  "  J  0Qsinu^t)    -    (y  Y5   "   5"  '^Qsina)2t)] 

+K4  |[(yi   -   Y^sintu^t)    -    (|y5   -  I  e^sinooj^t)    +    C^y^    -  -  ^^sini^2^)-\ 

~^1  H:(yi-YoSina3jLt)    +    (2   Y3-   ^  e^sinajj^t)    +    ( J  Ys"  -  $QsirD2^^^ 

L  L  L  W  W 

-^2  ^(Yi-^osinw-Lt)   +   (- Ya-  -  eQsinwj^t)    "    (^  Y5   "  y  *o^^^'^2^^^ 

roll  equation-- 

I2y6   =   -   KjL^(Yi-YQsina3^t)   +    (2  Y3-  J  eQsinoj^t)   +   (|  Y5-   j  $Qsinw2t)] 
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+K22-t(yi-YQ£into^t)    +    (j- 73"  ^  OQsinw^t)    -    (^75-    2    ^Qsinw2^)^ 

Driving  Force 

The    force   v.'hich   excites   the    cradle   is   transferred  through   the  springs 
to  the   cradle.      In   this   analysis   it   is   assigned  that   Lulu  is  heaving  and 
pitching   sinusoidal  with   the   same    frequency   W3_,    and  amplitudes   Yq   and   6^ , 
respectively.      Additionally,   it   is   assumed  that   Lulu  is   rolling  sin usoi dal- 
ly v/ith   amplitude   i     and   frequency   ^2. 

Drag  Forces 

Tnis   - 

low  Reynold's  numbers    involved  the   drag   forces    are   primarily   caused  by 
hydrodynamic   form  drag  rather  than   viscous    drag   and   are  nonlinear.      In   ad- 
dition,   for  geometric  shapes  with  hard  edges   the   drag   coefficients   are 
generally  independent   of  Reynold's  numbers. 

The  general  equation   describing   drag   forces   on  such   a  shape   is   given 
as     Drag  =     hp   A^  Cd  V    ,   and  it   is   in  this  where  we   find  the   problem.      If 
drag   is   to  alv;ays   oppose  motion   as    it   does  here  ,   the   equation    above   should 
actually  be  written      Drag   =  hp   A^Cd|v|v,  v;hich  will  preserve   the    correct 
sign. 
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Figure  6, 


J. n i v^ i s. J. j.y  8.n  £xs.Co  s.'^'CSuip'C  3.x.  uns  3wj.u.Gi0ii  wcs  iyi"j.Su..   j-ns-smiiCii  cc  3. 

differential  elerient  of  drag  is  defined  in  quadrant  (l)  of  Figure  6  as 

a(drag)  =  HpC-r)dy.dz(y+x9cos(|)COse+Z(i)co34)Cos9  )  |y+x9cos(J)COs9+z<})cos4)COse  |  . 

then  the  total  drag  on  auadrant  (l)  would  logically  "be 

D?AG  =  j  j   ^pCpjdxdz(y+x9cos9cos(J)+z6cos9cos(}))  |y+xecos9cos(})+Z(J)COs6cos(|)  j  . 
0  0 

This  approach  could  ce  pursued  to  an  exact  solution  if  an  identity  for 

(A+?+C/!  A+3+c!  '^as  obtainahle.   However,  (A+B+C)  |a+B+c|  only  equals 

(A+n+C)  (A+B+C)  when  (A+3+C)  is  a  positive  number.   The  determination  of 

:ho  e:-;a-t  sign  of  ( A+B+C )  or  (y+x9cos9cos(J)+z^cos6cos(J) )  can  only  be  obt.*ined 

vher.  a  >r.cwn  value  of  X,  or  Z  is  available.   Consequently,  if  X  and  Z  are 

ir.  *.-.e  ir.te.;;rand  they  cannot  be  assigned  a  specific  value  as  they  vary  from 
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0   to   L/2   and  V;/2   r..;i-,poc-'-i  vely .      Thus,    ar>   ir.aLho;natically   satic- lying   as    an 
exact   oquation    for   drag  v/ould  be ,  it    is   ala.s,  impossible   to  obtain.      Due   to 
this,   a  more   straight    for-^/ard  yet   more    tedious    finite   element   method  v/as 
used   to  approximate   drag   forces. 

The    cradle  was    first   divided   ud   into   four  Quadrants    as   shov/n  below: 


3  2 

h  1 


Figure    7. 

Using  the   sign   convention  previously   described,   the  velocities   of  any  point 
in   the   four  quadrants   are   as    follows: 


V      =    (y   +   X   9    cos9cos<J)   +    z({^cos8cos({)) 


V^    =    (y   +   X   9    cos9cos(J)    -    z(t)Cos9cos4') 


V      =    (y   -   X    9    cos9coS(|>   -    Z(i)cos9cos({)) 


^M   ~^y   ~   ^   9    cos9cos<^   +   z4)cos9cos(})) 


Each   quadrant   is   then   divided  into  sixteen   areas    and  the    drag   is    calculated 
for  each   one   as   shown   for  element  a       below. 
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Drag  element   SL22 

=  i^  pCj-j  Area   (a22)    Va22|Va22 
W        L 


V/ 

a 


Area  a 


22 


8       8 


"^22   ~    ^y^  ~T'  L9cos9coscJ)+  — r  V/(})cos6cos4) ) 


16 


16 


^11 

''^12 

ai3 

^-^ll. 

^21 

^22 

^23 

^21; 

^31 

^32 

^33 

a-3k 

^Ul 

a-U2 

^1^3 

^kh 

-^L/2        _    _ 

(Quadrant   l) 
Figure   8. 


3    . 


Drag(a22)   =  h  PCj^^o"  '    A'^^y*'  77L9cosecos(^+  ■^//4)Cos6cos(j))  ly+  7^0cos9cos(+) 

3  . 

+  -— ^/Jicosecoscbl 
16  ' 

Hence,  the  velocity  is  assumed  to  be  fixed  over  the  entire  element  and  thus, 
so  IS  "cne  sign  01  i,ne  arn.g.   xaeaxxy  ,  aa  une  mimutii-  ux  cxcmcxios  ^^^a    uj^ ,  Lii<= 
possibility  of  the  total  velocity  changing  sign  in  an  element  goes  down  un- 
til, in  the  limit,  you  have  numerically  integrated  the  vhole  surface  to  an 
exact  result.   In  this  analysis  6k   elements  are  taken  each  time  step  and 
added  together  to  produce  the  resultant  drag. 

Pitch  and  roll  drag  are  calculated  in  much  the  same  way  but  with  the 
added  complication  that  each  drag  element  acts  through  a  finite  moment  arm 
to  produce  the  drag  in  terms  of  foot  lbs.    One  difference  is  the  proper 
description  of  the  velocity  of  the  element  which  for  our  same  element  a22 
would  be 


Va22  =  (ycosScoscj)  +  ^  L*e+  -J-  V7^). 


3_ 

16 


The  computer  subroutines  "for  calculating  the  drag  forces  are  given  in 
the  Appendix. 
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In  cases  where  the  relative  velocity  between  a  submerged  body  and  its 
surrounding  water  does  not  remain  constant,  additional  forces  come  into  play 
as  a  result  of  the  acceleration  of  the  mass  of  water  entrained  by  the  body. 
This  added  mass  depends  on  the  physical  dim.ensions  of  the  body  as  well  as 
its  mode  of  motion.   The  actual  body  m.ass  plus  the  added  mass  is  called  the 
virtual,  induced,  or  total  mass.   Similarly,  there  are  added  moments  of 
inertia  terms  which  accompany  the  body's  angular  acceleration.   Conse- 
quently, for  a  long  flat  plate  of  dimensions  such  as 

I  1 


2a' 


the  added  mass  of  the  entrained  water  for  motion  perpendicular  to  the  plane 

W 
of  the  cradle  is  defined  as  Mg_  =  pira^,  or  in  our  case  a  =  '^i   aeptn  =  l,  so 


M  =  ^pttW^L.   The  total  mass  is  then  actual  mass  plus  added  mass 

Hp  =  M( actual)  +  M  (added) 

Mr  =  508  +  ^pTTW^L. 
The  added  roll  mor.ent  of  inertia  due  to  the  entrained  water  is  defined  as 
or  in  our  case         lap  -   l/Spira'^, 

Ia2  =  l/8pTr(|)'^L. 

The  total  roll  moment  of  inertia  is  then 

m^       1    W^ 
^2  =   I^  ^  8  ^^"16  L. 


Similarly,  for  pitch  it  is 


ML2    1    L*"  ^^ 

I    =  +  --  pTT  —7  V/. 

1    12    8    16 


These  forces  consisting  of  buoyancy  and  weight  are  constant  and  thus 
do  not  enter  into  the  equation  in  any  v;ay. 
Overview 

The  equations  developed  to  this  point  are  nonlinear  due  to  the  drag 
terms  tut  fortunately  can  be  solved  by  various  numerical  methods  of  inte- 
gration.  'The  equations,  arranged  in  state  variable  form  are  shown  in 
Figure  9.   It  is  important  to  note  that  the  drag  terms  shown  are  extremely 
complex  and  are  functions  of  five  of  the  six  state  variables.   More  specif- 
ically, the  heave,  pitch,  and  roll  drag  forces  are  functions  of  the  follovr- 
ing:  y(2),  vertical  velocity  of  cradle;  y(3),  the  cradle  pitch  angle;  y(^); 
the  cradle  pitch  angular  velocity;  y(5))  the  cradle  roll  angle;  and  y(6), 

dr>=.gs  are  different  functions  of  these  variables. 

The  forms  of  these  functions  have  been  explained  in  previous "para- 
graphs. 
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(K^*:-.g':<3->i;i.) 


-DRAG   (heave) 

f(y{2).y{3),y(M,y(5),y(6)) 

finite  elerient 


(KaL-KuL^KiL.K.L)  f  (,(2,  ,,(3)  .y^,  .y(  5)  ,,(6, ) 


(Kj.'tKi'J-iCW-KjH 


finite  elener.t 


-DRAG   (roll) 

r(y(2),y(3),y(i-),y(5),y(6)) 

finite  element 


(KiL+K2L-K3L-KI,L 


(KjL'+KtL'»KiL?»K2l ')  g 


(KiWL-KlVL-KjVL-Kj-.  L)  - 


(-KjW-K-'-'-KjHtK^U) 


(♦K3LW-Kl,LW*KiLW-K2HL) 


(Kj'KgtKj-fKI.) 


(KiW-ff:|.»-K;«-K3W) 


(K;L->KyL-"-<UyiiL)  (KiW-K2l'-K3V*Kl,V) 


(-K3L-KI.UKiI.*K;L)  (K3I,'*Ki,l'«KiL'>M')  (K3t«-K|,LH*yjWL-K;UL) 

PI.  ''I,  111. 


(Kim.-!!l,WL-It2in.+K3U  .) 


(KiV2»!t],«'»Kj-j'tK3W!) 


Figure  9, 
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chapter  iii. 

develop:-:£tit  of  liiifar  approxim.'^tiohs 

As  discussed  previously,  the  equations  as  developed  are  nonlinear. 
Unfortunately,  this  fact  results  in  the  necessity  of  a  "brute  force  ap- 
proach to  the  solution.   That  is,  the  well  developed  test  and  analysis 
methods  v;hich  apply  to  systems  of  linear  equations  cannot  be  used  to  draw 
■general  conclusions  about  a  nonlinear  system.   Consequently,  linear  ap- 
proximations for  the  drag  terms  would  be  helpful  only  if  they  were  shown 
to  give  reasonable  results  when  compared  to  the  exact  approach. 

Ideally,  the  nonlinear  drag  terms,  which  now  are  functions  of  five 
of  the  six  state  variables ,  should  be  reduced  to  functions  of  only  one 
velocity  term.   Actually,  it  would  be  reasonable  to  expect  the  majority 
oi    a  par uxuux'dr  urag  oexxii  i^u  "uc  uduacu.  uj     L,ln_  v^lw-^lL^  ^^l..l^  __-.-_--:l -'.-■' 
with  that  term.   That  is,  heave  drag  should  mainly  be  caused  by  heave 
velocity,  pitch  drag  by  pitch  angular  velocity,  and  similarly,  roll  drag 
by  roll  angular  velocity.  .    . 

Therefore,  using  this  as  an  initial  approach,  the  single  variable 
heave  drag  equation  might  be: 

Drag  (heave)  =  h     pCj)WLy2 
The  equation  for  pitch  drag  would  be  somewhat  more  complicated  and  would 
be  calculated  as  shown  in  Figure  10. 
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EiaLi 


IW^'^' 


Figure   10. 


d    (drag)   =   2  x  ^^  pCj5WdX(xe )  (xe  )X 


DRAG 


1/2 


4o2 


DRAG   (Ditch)    = 


pc^WL^e 
6h 


Similarly  the   equation  for  roll  drag   is 
DRAG    (roll)    = 


pCj)WL^(i>2 


b4 
As  yet,  even  these  equations  are  nonlinear  due  to  their  squared  velocity  terms. 

Several  approximations  to  velocity  squared  damping  have  been  used  in  the 
past.   One  of  the  commonly  used  electrical  approximations  consists  of  equat- 
ing the  energy  dissipated  by  a  current  squared  resistor  to  the  energy  dissi- 
pated by  a  linear  resistor  vhen  both  are  driven  by  a  sinusoidal  current 

source.  VThen  this  is  done  the  I^  term  is  in  effect  related  to  a  constant 

ft 
multiplied  by  a  single  power  of  I.   The  constant  which  results  is  —  -'-n^l 

where  I„  is  the  raaximuia  current  value  and  w,  is  the  driving  frequency. 

J?herefore,  if  we  consider  the  sea  as  a  velocity  squared  damper,  a  rea- 
sonable approximation  to  heave  drag  would  be        '     ' 


8 


J,pCi5WL(-  YQCo,)y 
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where  ^tYq-  is  the  rcaxiniun  velocity  of  y.   If  this  approach  is  applied  to 
pitch  and  roll  the  following  equations  result: 

h.OpC-QU-rfrpiy 


Drag  (heave)  = 


3tt 


Drag  (pitch)  =     2h-n 
Drag  (roll)   =■  ~ 


2UfT 

To  check  the  validity  of  this  linearization  approach  as  veil  as  the  ac- 
curacy of  the  6h   element  method,  a  comparison  vas  made  with  values  obtained 
from  the  exact  drag  equations  shown  in  Figure  10.   Since  these  exact  equa- 
tions are  only  valid  for  cases  where  single  velocity  terms  are  used,  the  com- 
parison was  done  for  a  case  where  heave  drag  is  only  a  function  of  heave 

drag  is  only  a  function  of  roll  angular  velocity.   In  other  words,  simple 
roll,  pitch,  and  heave  drag  were  tested  independently. 

The  results  of  the  comparison  shown  in  Figure  11  indicate  clearly  how 
the  finite  element  method  compares  almost  identically  with  the  exact  values 
as  computed  with  the  equations  of  Figure  10.   The  one  difference,  that  being 
that  the  linear  approximation  gave  lower  peak  values  as  they  appeared  more 
sinusoidal  in  shape  was  anticipated  because  we  are  equating  energy  dissipa- 
tion rather  than  peak  values . 

Additionally,  to  assure  that  the  one  term  approximation  will  give  the 
desired  drag  force  at  the  proper  time,  a  comparison  was  done  between  actual 
and  linear  drag  values  which  would  be  created  during  normal  cradle  oscilla- 
tions.  Figures  12  through  ih   compare  drag  values  with  Lulu  exhibiting  a 
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Figure  11. 
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ten  second  ?.0^   roll  and  pitch,  and  a  ^.25  ft.  heave  ampliturle.   The  corr;pari- 
son  shovs  that  for  this  particular  driving  frequency  and  amplitude  the  linear 
approxir.ation  is  in  phase  with  the  actual  values,  yet  the  amplitudes  do  not 
correspond.   Another  comparison  at  twice  the  frequency  of  the  first  confirmed 
that  vhile  a  single  variable  drag  term  approximates  the  timing  of  the  drag 
force  veil,  it  does  not  predict  the  amplitude. 

In  an  attempt  to  better  estimate  the  drag  forces  several  two  and  three 
term  approximations  were  investigated  but  the  results  shov/ed  that,  in  gener- 
al, complete  confusion  reigned.   Even  worse,  the  drag  forces  were  no  longer 
in  phase  with  the  actual  values. 

The  incorrect  drag  amplitudes  result  mainly  from  the  terms  in  the  coef- 
ficient of  proportion  of  the  linear  drag  term.   In  particular,  for  heave 

o 

C'^Z-o   tl.:.  -cr.rtant  ^  ^-'C--,'   dep'=»^''''=  v.^t-viiy  nn  w-.  .  the  freauency  of  the  cradle, 
and  Yq_  ,  the  maximum  displacement  of  the  cradle.   The  cradle  frequency  t^i 
can  be  assumed  to  be  the  same  as  the  driving  frequency,  but  Yq ,  the  maximum 
cradle  amplitude,  cannot  be  assumed  to  be  the  same  as  the  maximum  driving 
amplitude.   In  actuality,  Yq  of  the  cradle  is  the  maximum  displacement  of  the 
cradle  as  given  by  the  equations  of  motion.   However,  that  value  of  Yq  de- 
pends on  the  chosen  value  for  the  drag  term,  which  in  turn  is  a  function  of 
Yq.   Basically  we  have  an  iterative  problem  where  values  for  the  drag  terms 
must  be  initially  ass^jmed  and  comparisons  must  be  m.ade  with  the  actual  non- 
linear responses. 

To  provide  a  baseline  for  comparison  a  brute  force  analysis  was  con- 
ducted on  the  nonlinear  equations  to  determine  the  actual  response  or  gain 
of  the  system  at  various  input  frequencies.   Ten  separate  frequencies  were 
analyzed  and  the  steady  state  gains  were  plotted  versus  frequency.   The  same 


-36- 
analyais  .v/as  rer^oatod  but  this  tir;.e  u;:inp"  the  sinsJe  term  linear  drag  ap- 
proxi .'nation 3  in  p3.ace  of  the  accepted  6h   eleinent  drag  values.   This  analysis 
v;as  also  conducted  at  the  ten  different  driving  frequencies  vith  three  dif- 
ferent driving  ajnplitudes  . 

Figures  15  and  l6  show  clearly  the  differences  between  the  linear  and 
nonlinear  responses  at  various  driving  frequencies.   The  roll  response  curve 
is  ommitted  since  it  falls  just  short  of  being  on  top  of  the  pitch  curves. 
These  figures  show  that  the  linear  approximations  peak  at  the  same  frequen- 
cies as  do  the  actual  response  curves  but  the  gains  vary  widely  depending 
on  the  input  amplitude.   In  general,  the  linearized  response  curves  show  a 
condition  of  lesser  damping  than  is  actually  the  case.   However,  as  the  driv- 
'ing  amplitudes  go  up,  the  apparent  damping  increases  due  to  the  presence  of 

iVi£»      ~l  c-^~^--<      ■'n^       --^       J— -•-.-!_,_       T^J..,T,,        'v       ^        J_        J-V.- . -^-«>-'-- J-       —-"^      — ,„^-,^  „^->-  A   r^-^ 

one     ±dLr^&l       vcLJ.U5     ^±      U.i.xvj.rig,      cuiipJ.xuu.Ucr      Vxp|/       m      one      i_^i;;j- x  XCxCi^vy      '-'J-      ^^  w_i^^j.   w^uOij.  . 

-  Since  the  main  value  of  such  an  analysis  is  the  prediction  of  the 
resonance  peaks ,  and  since  the  linear  and  nonlinear  responses  are  almost 

W  CO 

identical  except  in  the  rangeO.5  Tj  to  2.0  —  ,  it  is  valuable  to  correct  the 

n        w^ 

linear  response  curves  by  increasing  the  apparent  damping  until  the  linear 
responses  con-form  to  the  nonlinear  responses.   This  can  be  done  by  multiply- 
ing the  previously  discussed  linearized  drag  terms  by  a  constant  depending 
on  the  magnitude  of  the  driving  amplitude. 

Comparing  the  linear  response  curves  to  the  simple  single  degree  of 
freedom  system  shown  in  Figure  IT,  we  find  that  the  gain  of  such  a  system 

at  resonance  is  given  by  the  equation  — -  where  —  is  the  damping  factor. 

Cc       Cc 

Using  this  equation  we  find  that  the  nonlinear  heave  resonance  peak  compares 
well  to  a  linear  system  with  a  damping  factor  of  0.193.  Similarly  the  pitch 
response  compares  well  with  a  damping  factor  of  O.I67.   After  finding  the 
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indicated  linear  values  of  damping  factors  for  the  linear  response  curves 
vith  various  input  amplitudes,  a  plot  of  drag  multipliers  versus  amplitude 
input  vas  made.   This  plot,  shown  in  Figure  l8,  indicates  the  factor  by 
which  the  drag  terms  must  be  multiplied  in  order  to  make  the  linear  response 
correspond  closely  to  the  nonlinear  response  for  each  input  amplitude.   For 
example,  for  our  normal  ^4.25  ft.  heave  and  10  roll  and  pitch  amplitude,  the 
heave  drag  term  would  need  to  be  multiplied  by  2.7  >  the  pitch  drag  tei-n  by 
^.2l^,  and  the  roll  drag  term  by  k.O. 

Summing  up,  linearization  of  the  nonlinear  drag  terms  is  a  complex  mat- 
ter.  As  in  all  problems  in  the  physical  sciences  it  is  never  feasible  to 
express  such  linearizations  with  absolute  certainty  and  with  the  inclusion 
■of  all  possible  variables.   In  fact,  to  treat  any  physical  problem  in  a 
laathcmatical  way  it  i'j  necessary  to  first  make  ascumpticnc  about  the  system'?, 
behavior  that  m.ay  not  be  entirely  correct.   In  this  case  drag  is  assumed  to 
be  linearly  proportional  to  velocity,  but,  as  is  later  discovered,,  it  is  also 
related  to  the  driving  amplitude  and  therefore  appropriate  correction  factors 
must  be  applied.   In  spite  of  the  difficulty  involved,  the  value  of  this 
linearization  is  to  provide  the  coefficients  for  the  linear  A  matrix  which 
can  then  be  operated  on  to  find  the  general  trends  of  the  system. 
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CHAPTER  IV. 
SOLUTION  OF  EQUATIONS 

The   complete  solution   of  differential  equations   may  be   split   into  two 
parts .   a  steady-state  solution   and  a  transient  solution  which    dies   out   at 
a  rate   depending  on   the   degree   of  damping,  mass,   and  time.      Our  main   interest 
here   is  with  the  steady-state  solution. 

Time   Domain   Analvsis 


The   tim.e   dom.ain   solution   of  the  equations   of  motion  was    found  by  the 
use   of  num.erical  integration.      The   integration  was   done   using  the  newly   in- 
stalled Interdata  computer  facility   at  M.I.T.'s   Mechanical  Engineering   De- 
partment.     The   same   department's   program  DYSYS ,    (DYnamic  SYstem  Simulation), 
was   used   Idl-    Lue    uluie   uuiualii   analvsis    of  tl".^   wcnlinccir   ccu-tior.c.      DY?'^^ 
solves   a  set   of  first   order  linear  or  nonlinear  differential  equations    and 
can  be  particularly   useful   for  obtaining   transient   as   well  as   steady-state 
solutions   to  nonlinear  time   varying  problems.      DYSYS  uses   a   fourth  order 
Runge-Kutta  integration   technique   to  obtain   an   approximate   solution  to  this 
equation-- 

Given  Y^ ,   the   value   of  the   Y   vector  at  some   value   of  t,   say  t    ,   then 

■  ..:..  0 

the  approximate  value  of  the  y  vector  at  t  +  At  is 


At 
Y  +  At  =  Yq  +  -^Kq  +  2K;l  +  2K2  +  K3) 


the  K's  being  equal  to 


^1   =   ^(^0   +   ^'   Yq   +  h      Kq) 
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K3   =    F(tQ    +   At,    Yq    +   K2) 

The  Runge-Kutta  !!ethod   is   derived  fi"'om  a  Taylor  Series   Expansion   about   the 
point   tQ    and  Y    .      Theoretically,   the   method  v;ill   converge   to   an   exact   solu- 
tion  as      t   approaches    zero...     The  method  has  been   shov;n   to  be   stable,   self- 
starting  and  to  have   a  built-in  error  check. 

Since   the  program  is   set  up   for  simultaneous    first   order   differential 
equations,   the   three  second-order  equations   must  be   expressed   as   six   first- 
order  equations.      This    can  be   done  by   letting 

y  =  yd) 

y(2)   =  yd)   =  y 

y(2)  =  y 
Similarly, 

yd)  =  y 

y(2)  =  y 

y(3)  =  e 

...  y(u)  =  e 

y(5)  =  (J) 

■  y(6)  =  (J), 

and  hence  y(l)   through  y(6)   are   the   six  state   variables  v/hich    completely 
define  the  system. 

The   computer  program  v;hich  was  written   to  effect   the   solution  is   given 
in  the   appendix.      The  program  has   as   subroutines   EQSYM,  which   indicates 
"Equation   Simulation''   as  well  as   three   separate   drag  subroutines  which  must 


be   called  each   time   step   and   calculates   the   drag  at   each   point    in  the 
program.      A  tirr.e   interval   of  0.2   seconds  v;as   selected   as    a  time   step   and 
gave   reasonable   results. 

Frequency   Domain   Analysis 

Since  not   only  specific  solutions    for  given   parameters  were   desired, 
a  frequency   domain   analysis   resulting   in   Bode   plots   was    undertaken.      Bode 
plots   clearly   illustrate   the   relative   stability   of  a  system,  with  gain   and 
phase  margins   often  being  defined  in   terms   of  Bode   plots. 

This   analysis  was    conducted   using   the   M.I.T.    Mechanical  Engineering 
Department   computer  program  MBODE.      This  program   calculates    frequency   do- 
main gain   and  phase  plots    for  linear  dynamic  systems    (with   constant   coef- 
ficients)  expressed  in  matrix   form  such   as 

X  =   Ax  +  Bu 

y  =   Cx  +  Du        "  .      • 

where 

A  =   the   system  A  matrix 

B  =  the   input  matrix 

C  =   the   output  matrix   (usually  diagonal  ones   or 

appropriate  weighting   factors) 
D  =   the   direct   coupled  output  matrix   (usually  zero). 

The  fact  that  MBODE  only  operates  on  linear  constant  coefficient  equa- 
tions is  the  reason  we  were  forced  to  develop  linear  approximations  for  the 
nonlinear  drag  terms   discussed  .in  section   III.      The   Bode   plots  which 


I 


resulted  from  this  method  of  analysis  shov;  clearly  the  proj;erties  of  the 
system  and  these  are  presented  in  section  V. 
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CHAPTER  V. 
Discussion  OF  RESULTS 

Time  Doniain  Response 

For  this  analysis  numerous  trials  were  run  at  various  damping  levels, 
input  frequencies,  input  amplitudes,  and  spring  constants.   The  majority 

of  the  trials  were  run  at  the  baseline  conditions  of  U.25  ft.  ship  heave, 

o  Ibf 

10  roll  and  pitch,  and  20  -—  spring  constant  with  a  10%  deviation  among 

the  four  springs.  The  10/J  deviation -is  the  method  by  which  the  nonideal 
nature  of  the  springs  was  modeled  and  is  the  vehicle  for  energy  transfer 
between  the  different  motion  modes . 

In  a  discussion  with  J.  Mavor  of  W. H.O.I,  it  was  learned  that  the  pre- 
viously mentioned  amplitudes  with  a  5-6  second  period  could  be  considered 

situation  of  heave,  pitch  and  roll  and  indicates  that  the  transmitted  ship 
heave  is  reduced  from  U.25  ft.  to  about  5  inches  while  the  roll  arid  pitch 

are  reduced  from  10°  to  2  .   A  total  effective  spring  constant  of  80-100 

Ibf 

is  required  to  accomplish  this  decoupling;  however,  the  analysis  il- 
lustrates the  fact  that  ship-cradle  motion  can  be  uncoupled  if  soft  springs 
are  used. 

In  Figure  20  the  pitch  and  roll  amplitude  have  been  doubled  and  the 
ship  has  a  slow  10-12  second  period  to  give  a  "worst  case"  example  of  pre- 
dicted cradle  response.   Although  this  decoupling  is  not  as  good  as  with 
the  shorter  period,  the  heave  has  still  been  reduced  to  less  than  one  foot 
while  the  roll  and  pitch  have  been  reduced  from  20  to  10°  and  7  respec- 
tively. 
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Frequency  Domain  Response 

Since  specific  response  curves  do  not  reveal  much  about  the  general 
nature  of  the  system,  a  frequency  domain  analysis  was  conducted  in  order  to 
get  more  descriptive  information  on  how  the  system  will  react  to  various 
input  frequencies.   Figure  21  is  a  plot  of  the  system's  gain  response  at 
baseline  conditions  of  amplitude  and  spring  constant  hut  at  ten  different 
input  frequencies.   It  is,  in  effect,  a  Bode  plot  of  the  nonlinear  case. 
The  plot  shows  clearly  the  three  reasonances  of  the  system,  the  leftmost 
peak  being  the  heave  resonance  and  the  rightmost  peaks  being  pitch  and  roll 
resonances,  the  roll  resonance  being  the  smaller  of  the  two. 

This  figure  required  ten  separate  solutions  of  the  nonlinear  equa- 
tions of  motion  involving  great  expenditures  of  computer  time.   Although 
the  ulot  shows  the  response  of  the  system  under  fixed  condixions  of  spring 
constant,  mass,  and  driving  amplitude,  it  is  clearly  not  a  desirable  way  in 
which  to  find  how  the  system  responds  when  differing  sea  conditions  exist 
or  when  the  system  parameters  are  modified.   It  is  here  that  the  value  of 
the  linearization  process  comes  into  full  view. 

In  Chapter  III  it  was  noted  that  the  nonlinear  drag  terms  could  be 

linearized  once  the  driving  amplitude  was  known.   Again  using  the  U.25  ft. 

o 
heave  with  10  roll  and  pitch  as  the  baseline  we  can  find  the  linear  drag 

approximation  by  referring  to  Figure  l8.  Chapter  III,  and  obtaining  the 

proper  multiplier.   Once  these  drag  coefficients  are  set,  Bode  plots  may 

be  initiated  which  show  the  response  of  all  outputs  versus  all  inputs  over 

a  wide  frequency  range.  Additionally,  the  effects  of  variations  in  spring 

constant  or  mass  may  be  studied  with  only  minor  modifications  to  the  now 

linearized  A  matrix.  •. 
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Fig^^ss  22  throuG^i  25  show  the  system's  response  across  a  range  of 
.01  to  1.0  Hertz  or  a  roll  pitch  and  heave  period  of  1.0  to  100.0  seconds. 
The  outputs  for  heave  versus  ship  pitch  and  roll  show  that  the  ship's 
pitching  and  rolling  couples  energy  into  heave  motion  which  then  peaks  at 
heave  natural  frequency  as  well  as  the  pitch  and  roll  natural  frequency. 
This  fact  indicates  the  dangers  of  designing  a  decoupling  system  to  de- 
couple only  heave.   It  is  clear  from  the  plot  that  a  decoupling  system  de- 
signed only  for  heave  may  well  excite  the  pitch  and  roll  modes  at  their 
natural  frequencies.   This  fact  was  overlooked  "by  single  degree  of  freedom 
analysis. 

Several  other  conclusions  may  be  drawn  from  Figures  22  through  25. 
While  a  20  Ibf/ft  was  chosen  as  a  baseline  for  the  spring  constants  and 
was  Used  for  Figures  22  and  23,  it  v.'ould  be  desirable  to  predict  easily 
what  would  happen  with  higher  or  lower  spring  constants .   Figures  2U  and 
25  show  the  response  of  the  system  with  a  spring  constant  of  four  -times 
that  of  Figures  22  and  23.   It  is  not  surprising  to  see  that  the  resonant 
frequencies  are  shifted  upward  by  a  factor  of  approximately  two.   This  in- 
dicates that  in  spite  of  the  nonlinearity  of  the  parent  system,  the  system 
obeys  the  basic  laws  for  calculating  natural  frequencies  of  translational 

and  rotational  systems.   Being  more  specific,  heave  natural  frequency  is 

/Keff 
proportional  to    „    ,  while  pitch  and  roll  natural  frequencies  are 

proportional  to  v4x"  and    /^       respectively. 
II         12 

As  a  result  of  the  above,  we  have  information  on  the  system's  behavior 

over  a  wide  range  of  frequencies  without  the  need  of  calculating  a  Bode  plot 

for  each  value  of  K.   In  effect,  using  K  =  20  Ibf/ft  as  the  starting  point. 


-56- 
we  need  only  shift  the  response  peaks  up  or  down  the  frequency  scale  by 


'^ew        , 

a  factor  of  the  square  root  of  K^ev/^'20'  °^   w ^  "^"^old  ^new   • 

This  same  method  of  peak  shifting  can  he  applied  to  variations  of  mass  or 
moment  of  inertia  in  a  like  manner,  but  in  this  case  the  shift  will  be 

,    ^      /Mtlbasej  /l(base) 

proportional  to     7777 r-         or         — ;^ —     . 

^  ^  iMtCnew)  I(new) 

Sxwiming  up,  this  analysis  has  shown  that  motions  of  a  submerged  cradle 
can,  in  fact,  be  nearly  decoupled  from  surface  ship  motions  if  a  soft  ef- 
fective spring  constant  is  provided  at  the  four  cradle  corners .   In  addi- 
tion, the  analysis  has  shown  that  the  critical  design  consideration  is  the 
pitch  and  roll  natural  frequency  rather  than  the  natural  frequency  of 
hc:i"e.   '^c  °  vo=;irit. .  rareful  thought  must  be  given  Lo  the  selection  cf 
spring  constants  since  the  equations  of  motion  are  coupled  and  are  capable 
of  transferring  energy  from  one  mode  to  another  and  hence  exciting  other 
modes  at  their  natural  frequencies . 


I 


-57- 

aiAPTER   VI. 

APPLICATION   POSSIBILITIES 

It   is  not  the   major  objective   of  this   thesis   to  offer  a   detailed   de- 
sign of  a  decoupling  system.      It   is   enlightening,  however,   to   discuss 
several  approaches   that   could  be   applied  in   the   design   of  such   a  system. 
The  primary   consideration  here   is   to  obtain   a  soft  spring   constant   and 
still  support   a  large   static  load. 

From  the  previous    analysis   it   is   apparent  that  the  soft  spring   con- 
stant must  be   at   the    four  cradle    comers.      Further  review   indicates   a 
spring  constant  bet^^/een   20   and  50    Ibrr/ft  would  seem  to  be   a  reasonable 
selection   for  these   comer  springs. 

At  this   r^cir.t  th'^'>''°   ^'•^   f^o  rvDes   of  approaches   to   cout>ider.      The 
first  and  simplest   approach   is   the  passive   system  where  we   try   to  create 
a  spring  with   the  proper  softness   so  that  without   any   outside   control  the 
springs  will  decouple   the   ship-cradle  motions.      This    approach   can  be 
likened  to  vibration   isolation   of  a  motor  or  a   vibrating  mass.      It   is   sim- 
ple and  reasonably  foolproof. 

The   alternative  method  is   of  course   a  more    complex  active   system 
where   some   sort   of  control  is   used  to  actively  uncouple   the   two  motions. 
This   control   can  be   in   the   form  of  position,   velocity,   acceleration,    inte- 
gral or  differential  feedback,  horizontal  sensing  transducers   or  mechani- 
cal linkages.      Consequently,   unlike  passive   systems,    these   types   of  sys- 
tems  can  be   tuned  to  changing  sea   conditions   or  conditions   of  loading, 
thereby  giving  them  a  possible   advantage. 
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Passive  Syntoms 

To   create   a  passive   system   for  the    cradle   v;ould   require    consideration 
of  several  ir.ethods   of  achieving   the   desired   soft   spring  constant.      A   first 
choice    for  component  selection  would  probably  be   simple  metallic  springs 
v;hich   deform  in   accordance  with  Hooke's    Law.      Unfortunately,   the   static 
cradle  v;eight   of  about    16  ,000   pounds  would  stretch   the   four  50   Ib^/ft 
springs   to  about   80    feet  before   reaching  equilibrium.      Clearly,   this  would 
be   unacceptable   for  if  coil  springs  v;ere   \ised  they  would  necessarily   re- 
semble  long   "slinkys"'   and  would  prove    completely   unmanageable.      Leaf  springs 
would  similarly  be   of  little   use   as   they  too  would  deflect  excessively  in 
the  static  mode. 

Torsional  springs  would  possibly  work   if  they  were   rigged  through   a 
com.plex  geai'ing  system  v/hich  wound   up   as    load  was    applied;  however,   they 
too  would  require   a  large  static  deflection   under  the   simple   cradle   load 
and  would  involve   complex  machinery   and  cable   drums.    Such   a  mechanical  sys- 
tem would  probably  be  more   of  a  liability  than   an   asset   in   an   ocean  environ- 
ment. 

Another  possible   alternative  might  be  to  take   advantage   of  the  non- 
linear stress   strain  properties   of  elastomers.      Vulcanized  rubber  will 
stretch  to  nearly  ten  tines    its   own  length   and   support  up  to  ten  tons   per 
square   inch.      Unfortunately,   the  resistance  to  deformation  usually  increases 
in  greater  proportion  than  does   the   extension.      If,  by  special  treatment,   an 
elastomer  with  the   opposite  properties  were   developed,   it  might  be   a  logical 
choice  to  use   as   a  spring.      However,   a  short   literature   search  turned  up  no 
suggestion  that   such   a  material  exists. 


-59- 

A  more  rear.on'ible  method  of  producing  such  an  effective  spring  constant 

would  "be  vith  a  hydraulic  or  pneumatic  device.   A  device  such  as  a  hydraulic 
or  pneamatic  cylinder  has  the  advantage  that  it  can  be  biased  to  an  initial 
static  load  and  then  oscillate  about  that  equilibrium  position  with  an  ef- 
fective spring  constant  which  depends  on  pressure,  temperature,  properties 
of  the  fluid,  and  dimensions  of  the  system. 

To  illustrate  such  a  system  consider  Figure  26. 
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Figure   26. 
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This  system  would  have  a  spring  constant  which  would  depend  heavily  on  k, 
the  compressibility  of  the  fluid.   It  is  evident  that  a  fluid  with  a  high 
compressibility,  such  as  air,  would  result  in  a  softer  overall  spring  con- 
stant which  is  desirable  in  this  case.   Additionally,  a  decoupling  system 
with  air  springs  v;ould  be  capable  of  being  set  to  an  initial  static  force 
equal  to  that  of  the  cradle  weight.   This  might  be  accomplished  with  a  sim- 
ple pressure  regulating  device.   Following  initial  set,  the  cradle  could 
then  oscillate  about  equilibrium  with  a  soft  effective  spring  constant  de- 
pending on  the  previously  mentioned  parameters. 
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Beforo  considering  this  simple  approach  there  are  two  important  con- 
siderations to  be  a.vare  of.   First,  the  oscillating  piston  must  look  exact- 
ly like  a  soft  spring  to  the  disturbing  ship  motion.   That  is,  the  differ- 
ence in  piston  force  from  eq.uilibrium  to  the  extremes  of  piston  travel  must 
compare  well  with  what  would  be  seen  if  a  real  metallic  spring  were  being 
used.   In  other  words,  for  a  ten  foot  pneumatic  cylinder  with  a  normally 
centered  piston  (capable  of  ±  5  foot  excursions),  the  fully  extended  piston 
force  must  be  equal  to  static  cradle  weight  plus  K  (effective)  "  5  feet. 
For  example,  if  the  desired  spring  constant  were  50  lb|>/ft,  then  the  piston 
could  only  exert  250  Ib^  more  than  the  cradle  weight  at  its  maximum  stroke 
of  five  feet.   Quickly  calculating  the  change  in  pressure  for  full  compres- 
sion using  P]_V2_  =  ^2'^o-,   i't  is  obvious  that  if  the  pneumatic  system  were 
dead  ended,  the  reduction  in  voluiae  Lo  1/5  of  its  previous  value  would  re- 
sult in  a  fivefold  increase  in  pressure. 

Consequently,  the  first  consideration  in  such  a  design  is  the  fact  that 
a  plenum  chamber  would  have  to  be  used  in  conjunction  with  such  a  system  in 
order  to  have  the  desired  results.   To  illustrate  the  size  of  the  required 
plenum  chamber,  we  can  again  use  P^^-i  =  P2^2  ^^^  "^1  ~  Tg).   Assuming  a  base- 
line piston  area  of  1  ft.^,  a  static  cradle  load  of  i+000  Ib^  per  spring 
would  require  a  piston  pressure  difference  of  278  psi  in  order  to  maintain 
the  cradle  at  equilibrium.   Calculating  the  plenum  size  required  with  the 
following  conditions 

AV  at  full  stroke  =  5  ft. 3 
AF  at  full  stroke  =  250  lb» 
^1^1  =  ^2^2    o^     278  V-L  =  (278  +  1.7)  (V^  -  5  ft^) 

V3_  =  815  ft^ 
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ve  unfortunately  find  it  is  a  chamber  far  exceeding  the  capacity  of  the 
cradle  to  easily  carry. 

Hov/ever ,  further  calculations  show  that  the  size  of  the  required  plenum 
will  go  down  as  working  pressure  goes  up.   For  example, 
assume     P  =  2000  psi 

Piston  area  =  2  in . ^ 


or 


2000  V^  =  (2000  +  125)  (V^  -  .O69) 


and  is  a  volume  which  is  well  within  the  carrying  capacity  of  the  cradle. 

The  second  major  consideration  in  passive  system  design  is  the  problem 
of  centering  the  system.   With  a  pneumatic  spring  such  as  is  sketched  in 
Figure  27,  the  piston  and  pleniom  might  be  pressurized  on  the  surface  until 
the  piston  was  positioned  at  midstroke  and  supporting  the  cTouile's  weight. 


PISTON 


Figure  2? 
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Tiiis  initial  setting,  although  correct  on  the  surface,  would  be  incorrect 
after  the  effects  of  :;ydrostatic  pressure,  temperature,  and  leakage  had  taken 
place.   This  fact  vould  necessitate  the  constant  bleeding  in  of  make-up  gas. 
In  effect,  sone  active  control  would  be  necessary  to  make  this  passive  design 
workable. 

In  addition  to  the  aforementioned  major  design  considerations,  when  us- 
ing a  gas,  with  its  higher  compressibility,  there  are  serious  technical  prob- 
lems to  be  overcome.   These  problems  involve  containment  of  a  gas  under  high 
pressure,  sealing,  and  danger  of  explosion ,   all  of  which  may  lead  to  bulky 
and  expensive  design. 

Throughout  this  entire  first  section  only  passive  decoupling  systems  are 
discussed.   These  systems  all  have  the  same  characteristics  in  that  they  os- 
r.Tnnt.o  i^hm^t   an  t=-aun  ibri-am  "Dositicn  not  determined  actively  by  sensors  ejid 
feedback,  but  by  cradle  weight,  system  dimensions,  and  initial  conditions. 

Active  Systems 

Systems  which  employ  active  feedback  for  control  are  called  active  sys- 
tems .   Typical  examples  of  these  systems  include  automobile  power  steering , 
hydraulic  positioning  devices,  and  automated  milling  machine  tables.   These 
active  systems  are  all  characterized  by  a  mechanical,  electrical,  or  hy- 
draulic system  by  which  information  on  the  output  is  fed  back  to  control  the 
actions  of  the  input.   This  feedback  control  can  be  designed  to  provide  fast 
or  slow  response,  maximum  or  minimum  overshoot  or  undershoot,  and  nearly  any 
natural  frequency  within  the  constraints  imposed  by  the  parameters  of  the 
system. 

The  simplest  form  of  an  electrohydraulic  position  control  device  is 
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shown  below   in   Figure  28. 
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Figure  28. 

The  reeabacK  loop  consists  oi  sensors  whicn  leea  uacis.  pooxoxoii,  vcj.^-il— 
ity,  and  acceleration  information  which  improve  the  system's  response  to 
disturbances  at  the  input.   Position  feedback  often  is  sufficient  for  system 
control  but  does  not  provide  the  tight  feedback  loop  that  results  when 
velocity  and  acceleration  are  also  used. 

A  typical  response  curve  for  such  a  system  is  shown  in  Figure  29- 
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That  is,  the  Gystcm  vrould  tend  to  rise  to  the  desired  level,  overshoot  it, 
and  then  oscillo.te  about  that  level  with  a  decaying  oscillation  depending 
on  the  para-Tieters  of  the  system.   It  is  this  oscillation  which  might  well  be 
used  to  create  a  virtual  soft  spring  constant  at  the  constant  load. 


A  Proportional  Control  Approach 

As  previously  stated,  feeding  back  all  three  variables  would  tend  to 
make  the  system  have  less  overshoot  and  oscillation.   Since  we  are  attempt- 
ing to  utilize  the  system's  natural  oscillations,  we  would  be  more  success- 
ful using  only  a  position  feedback  system  which,  while  not  enhancing  oscil- 
lations, does  not  discourage  them  strongly. 

To  determine  the  feasibility  of  such  a  system  and  to  get  some  design 
parameters  for  guidance,  the  system  equations  are  determined  below. 

First,  looking  at  the  spool  valve  and  assuming  the  hydraulic  fluid  is 
water. 

Valve  flow — 

•^m  =  ^ox  ^  "  ^op^m 
where  Q^^  =  flow  to  load 

Kqx  -  flow  sensitivity  or  valve  gain 

K       =  negative  slope  of  valve  characteristic  curve  at 

constant  value  of  x.   See  Figure  30. 


Figure  30. 
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For  an  ideal  valve  the  slope,  K   ,  is  very  snail  and  can  be  assumed  as  zero. 
Such  a  valve  is  called  a  flow  control  valve. 
Valve  forces — 
-KrX  +  Fj^  =  0 

where  Kp  =  valve  control  spring 
F.  =  valve  force 
Next  5  considering  the  piston, 
Piston  flow — 

^  =  '•*  -^  Pm  '-5  (Ke  -H  Ij   *   CgP^ 
where : 

C2  =  piston  leakage 

Kg  =  compliance  of  piston 

V  =  volume  of  piston  chamber 

Bjjj  =  bulk  modulus  of  fluid 


Piston  forces — 

M"'  =  P  A  +  \L 
y    m     t 

A  =  area  of  piston 
M  =  mass  of  load 


W  =  weight  of  static  load 


Combining  all  four  equations, 

m 

and :   ^!y  =  P  A  +  W+  . 
**    m     t 

Finally, 


-jihiK,^-))   ^-^.  Ay  =  K^^X 


"m 
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and  using   Lr^Place  tranc forms, 

K^,,X(s) 


y(s)  = 


^-•(V  4-   _1   4-   c   ^ +   ^j 


s(s2-^(K  +  —-)  +  s 

2A      ^m 


A 


Such  a  system  can  "be  described  in  terms  of  poles  and  zeros  and  valuable 
information  can  be  obtained  with  an  S  plane  plot  of  the  position  of  these 
poles  and  zeros.   Our  case  results  in  three  poles,  one  at  the  origin  and  a 
pair  of  complex  conjugate  poles.   The  conjugate  pole  pair  indicates  the 
system  is  oscillatory  while  the  pole  at  the  origin  indicates  that  the  system 
may  not  reach  the  desired  new  setpoint  or  that  it  will  have  a  steady  state 
error.   However,  this  is  of  little  concern  to  us  here.   Plotting  the  zeros 
in  Figure  31,  ^y^i 
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Figure  31. 


The  natural  frequency  of  interest  is  defined  by  the  length  of  the 
phazor  connecting  the  origin  to  a  complex  pole.   Basically,  if  we  put  the 
transfer  function  in  the  form 


T(s)  = 


const 


s(s  +  a^  +  JS-l)(s  +  a-L  -  jB-j_) 
we  can  adjust  the  natural  frequency  of  the  system  by  adjusting  the  coef- 
ficients which  control  the  position  of  the  complex  conjugate  poles. 


I 
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In  our  case, 


M      V 
calling   ^A^f^e  +  T^  =  ^2 


M 


A  ^  1 


=  ^0 


and  revritins 


y(s)  = 


const 


s(A2s2  +  A-^s  +  Aq) 


The  natural  frequency  of  the  system  is  defined  by  the  roots  of  the 
equation  ,2  .  :^  „  .  M 


or  the  magnitude  of 


or 


s"^  +  7—  s  +  —  =  0, 
A2     A2 


At  _  /,An  ,2  ,  ,_ .  ,An, 

A2  A2  A2 

/  Ai  2     At  2     aJJ 


If  we  assume  leakage  is  equal  to  zero  then  Cp  =  0  and  all  terms  with  A 
go  to  zero,  and  the  expression  reduces  to 


~=Wa 


0)   = 


M^(Ke  +  -)) 


=/ 


2A2 


e    3 

To  get  a  natural  frequency  of  0.03HZ  or  ahout  0.2  the  argument  of  the 
square  root  would  have  to  be  about  O.5O. 

Assaying  a  baseline  piston  size  of  1  ft^  and  V  =  10  ft^  using  a  mass  of 
2000  slugs  and  the  bulk  modulus  of  water  3=  320,000  psi,  the  equation  gives 
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2  ft^ 


^'^^   sec^     2000  slugs  {K^   +    ^^^    )   , 


1_ 

^  ~   2000  slugs  (K^  + 

32,000 

or,  vith  proper  units 

1  (im4)(liiM(2)    in^ 


0.50 


3eo2    -     ,000  l^f  -   ^^^1       _ft_(K  ,   1T280   in3     ^ 
ft  12in     ^      320,000psi 


•     3 

K„  ^   1497.00     ^       . 
s  psi 

Hence  a  compliance  of  about  500  rrr"  must  be  built  into  the  system.   Assuming 
an  air  accumulator  with  a  working  pressure  of  278  psi  and  change  of  one  psi 

P^V-L  =  P2V2 

2T8V^  =  279  (Vj^  -  500) 

V^  =  80  ft^  , 

•i.tVit  r*!-!     oo-a-in     -t  c;     ve^r-v     "nrcp     pnn     "nnlVv. 

Fortunately,  as  was  shown  in  the  passive  systems,  this  large  accumiUlator 
may  be  made  smaller  by  use  of  a  higher  working  pressure.  For  example, 
assume      P-,  =  2000  psi 
Area  =  2  in^ 

2(li)  in^ 
0.50  =  


2000  slugs  1ft  ,  20 

12  in      ^  e   320,000  ) 

3 


K  -  0.096  ^ 
^  psi 

or  again,  using  an  air  accumulator  at  2000  psi,  one  psi  should  change  the 

volume  .  3 

0.096  '" 


psi 
PlVl  =  P2V2 
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2000  V^  =  2001  (V^  -  0.096) 


192  in^  =  V^ 

which  is  of  reasonable  size  and  is  easy  to  build. 

Actually,  a  working  pressure  somewhere  between  the  limits  discussed 
here  would  be  most  appropriate  in  order  to  reduce  piston  size  and  weight 
as  well  as  reduce  the  plenum  size.   In  short,  there  is  probably  an  optimum 
selection  of  working  pressure,  piston  area,  and  hydraulic  fluid.   A  no- 
tional system  employing  this  approach  is  shown  below  in  Figure  32. 
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Figure  32. 
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A  Proporticnal-Into^ral  Control  Approach 

By  increasing  the  level  of  control  is  include  integral  as  veil  as 
proportional  control  we  may  be  able  to  improve  the  system  and  possibly  re- 
duce the  size  of  the  plenum  chamber. 

Modeling  the  system  as  shown  belov;  the  following  equations  result. 
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Figure   33. 
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Combining  equations  (l),  (2),  and  (3) 


1^  =  -  P  -  Fovsinwt  -  Dy 


but 


•    _y_  ^  Qm 


■m   Ket   AKet    * 


Plugging  in 


Area  =  A 

Pj^  =  Load  pressure 

^et  -   Total  fluid  plus 

cylinder  compliance 

Qj^  =  Flow  to  load 


FQsinwt  =  Disturbance 


trol  constant 


S  =  Integral  control 
constant 

D  =  Damping  constant 


Qm 

(1)  y  =  K^t^m  -  1 

(2)  My  =  -Pj^A-FoSinwt-Dy 

(3)  Qm  =  Ry+s/ydt-c(t). 


my  =  -  k{- —  +  Tv~r)    -   Fowcoswt 
^et   -^^et 


-  Dy 
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al2o, 


Q^_  =  Ry  +  Sy  -  C(t) 


or  finally,  p^     ^^  ^      ^  Sy 

•  l^'    +  Dy  +  Z — y  +  - — y  +  =  C(t  )+Fov2sina)t . 


Using  LaPlace  transforr.s, 


_  C(t)    +  Fow^sincot 

y(s)   = 


Ket  Ke^       K^t 

Where  C(t)    is   the   input   setpoint   (time   rate   of  change)    and  Fov^sincot   is  the 
disturhance ,   the   characteristic   equation   is   then 

^et      ^et    et 

or  T>o-:i      „         ^^       „ 

S^  +  — —  +  7,-^  b-  + 


M    KetM      Ket^^   KgtM 
This  equation  can  then  be  analyzed  with  standard  procedures  and  proper  co- 
efficients for  R,  S,  and  Kg^  can  he  found. 

Some  general  comments  can  be  made  about  proportional-integral  versus 
simple  proportional  control.   Addition  of  integral  control  allows  the  sys- 
tem to  respond  to  a  change  in  set  point  with  zero  steady  state  error  while 
it  does  not  change  the  basic  natural  frequency  of  the  system.   Moreover, 
this  control  element  provides  an  additional  independent  constant  which  can 
be  used  to  more  accurately  tune  the  system  for  maximum  decoupling  effect. 
The  availability  of  this  additional  constant  may  well  offset  the  complexity 
at^ded  with  an  integral  feedback  unit . 

One  fact  to  be  aware  of  with  the  P.I.  control  discussed  above  is  that 
the  cradle  has  been  shown  to  exhibit  ].ess  than  critical  damping.   Hence, 
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the  chara.cterli:tic  equation  would  probably  become  unstable  upon  the  ap- 
plication of  a  reasonable  gain  and  additional  damping  may  be  needed.   ITiis 
damping  could  be  provided  either  by  mechanical  means  or  by  additional  feed- 
back such  as  velocity  and  acceleration. 

The  determination  of  the  optimal  design  for  such  a  feedback  control 
system  is  left  for  future  work. 

Other  Systems 

Basically,  all  other  systems  are  combinations  of  the  passive  and  active 
systems.   Some  system  modifications  which  may  be  desirable  include: 

(1)  Inclusion  of  horizontal  sensing  valves  which  tend  to  channel  fluid 
to  the  low  sides  of  the  cradle  thus  giving  horizontal  feedback  to  an  active 
system. 

(2)  Substitution  of  rotary  air  or  fluid  motors  vith  gear  racKs  to 
provide  the  vertical  position  control. 

(3)  The  utilization  of  an  active-passive  system  which  employs  a  posi- 
tion control  to  center  a  passive  piston  or  rack. 
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CHAPTER  VII. 

C0IICLUSIC::3  mid  RECC:"!E.IDA?IOriS 

As  a  result  of  the  investigations  carried  out  in  this  thesis,  several 
conclusions  have  been  arrived  at. 

1.  A  four  point  passive  decoupling  scheme  can  substantially  decrease 
ship  heave  5  pitch,  and  roll  transmission  to  a  submerged  cradle. 

2.  The  effective  spring  constants  required  for  decoupling  fall  in  the 
range  of  20-50  Ib^/ft. 

3.  Each  mode  of  oscillation  will  couple  energy  into  the  other  tvo 
modes,  hence  the  design  of  a  decoupling  system  must  take  each  natural  fre- 
quency into  consideration. 

k.  Cradle  hydrodynamic  drag  is  a  complex  function  of  frequency,  ampli- 
tude, and  velocity  ana  musi-  De  cnosen  x,o  suxu  une  j.iiux^xu.ua.1  olouicxolv^x^  .,1.1  ;li 
is  to  be  modeled. 

5.  Conventional  springs  such  as  coil  or  leaf  springs  which  will  have  a 
soft  spring  constant  and  also  a  large  initial  static  load  do  not  exist. 

6.  Passive  pneumatic  springs,  while  meeting  desired  requirements,  would 
require  large  plenum  chambers  in  order  to  operate  at  safe  and  reasonable 
pressure  levels.   Additionally,  these  passive  devices  would  be  difficult  to 
center  in  the  cylinder. 

7.  A  logical  approach  for  creation  of  soft  springs  with  initial  static 
loads  would  be  an  active  feedback  control  system.   A  system,  such  as  one 
utilizing  proportional  and  integral  feedback,  could  substantially  decouple 
ship  motions  from  the  cradle,  yet  maintain  the  resultant  piston  oscillations 
about  the  center  of  the  cylinder  or  set  point. 
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In  summary,  both  the  theoretical  and  practical  aspects 
of  a  cradle  mounted  decoupling  system  have  been  investigated. 
The  value  of  this  analysis  is  that  it  brings  us  to  the  point 
where  we  must  decide  whether  the  possible  improvements  in  op- 
erating conditions  warrant  the  considerable  time  and  expense 
involved  in  making  such  a  system  come  true.   If  the  analysis 
had  shown  that  a  simple  passive  system  V70uld  work  and  was 
readily  available,  further  tests  and  construction  would  be 
easy  to  justify.   On  the  other  hand,  since  an  active  system 
seems  to  be  necessary,  extreme  caution  should  be  exercised 
prior  to  construction  to  assure  that;  first,  proper  feedback 
coefficients  are  found;  and  second,  the  necessary  control  com- 
ponents are  available  r.nd  are  reaaiiy  aaapraoie  -co  cue  octictii 
environment.   In  addition  to  the  above,  care  must  be  taken  to 
insure  that  any  resultant  detailed  design  satisfys  all  safe- 
ty requirements  which  may  be  imposed  by  WHOI,  the  Navy,  or 
any  other  regulatory  bodys . 

Other  means  of  decoupling  ship-cradle  motion  have  been, 
and  are  being  evaluated.   The  author  hopes  that  this  work  will 
enable  WHOI  to  more  clearly  evaluate  the  cradle  mounted  system 
and  how  it  compares  to  other  methods  being  considered. 
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